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Introduction 

Fiber-reinforced  composites  are  finding  use  in  increasing  numbers  of  aerospace 
applications.  One  reason  for  this  trend  is  that,  from  a  structural  perspective,  these 
composites  offer  outstanding  strength-to-weight  advantages  over  traditional  materials 
such  as  metal  alloys.  A  common  problem  in  composites,  however,  is  the  inability  of 
current  methods  to  accurately  predict  their  strength  and  durability.  The  innate 
heterogeneity  on  the  constituent  level  creates  a  complicated  strain  field  and  can  create 
stress  concentrations  in  the  composite.  The  type  of  composites  studied  within  is  fiber 
reinforced  composites,  but  the  formulations  extend  to  fully  three-dimensional  composites 
such  as  woven  composite  architectures.  Because  of  the  size  and  shear  number  of  fibers  in 
these  composites,  it  is  not  practical  or  computationally  possible  to  explicitly  model  the 
fiber  and  matrix  material.  Thus,  homogenization  methods  are  often  used  to  replace  the 
heterogeneous  composite  with  a  homogeneous  block  of  material  having  representative 
mechanical  and  thermal  properties.  By  homogenizing,  however,  the  local  details  of  the 
microstructure  are  lost.  It  is  these  local  details  that  can  play  large  roles  in  stress 
concentrations  and  ultimately  failure  of  the  composite.  De-homogenization  techniques 
provide  the  recovery  of  these  important  local  details. 

One  type  of  stress  concentration  occurs  when  the  microstructure  changes  abruptly 
(e.g.  a  ply  interface).  The  first  goal  of  this  work  was  to  create  an  analytical  tool  to  model 
and  predict  the  strain  felt  on  the  level  of  the  microstructure  in  the  composite  where  the 
microstructure  abruptly  changes. 

Task  1:  Formulation  of  the  asymptotic  theory  for  abruptly  changing 
microstructures 

It  was  shown  in  [1]  that  periodic  homogenization  methods  are  sufficient  to  predict 
strain  states  when  the  microstructure  is  uniformly  periodic.  Figure  1  illustrates  the 
success  of  these  methods  on  a  free  edge  problem  comprised  of  a  0/90/0  IM7/5250-4 
composite  laminate  subjected  to  a  1%  displacement  strain  in  the  0°  direction.  The  failure 
criteria  J1  and  J2  are  dilatational  and  deviatoric  strain  invariants  known  in  the  literature 
as  Strain  Invariant  Failure  Theory  (SIFT),  introduced  by  Gosse  et  al  in  [2],  The  curves 
referred  to  as  “Ply  RVE”  in  Figure  1  are  the  curves  computed  according  to  the  periodic 
homogenization  in  [1]  and  the  curves  referred  to  as  “3D  Sim.”  are  the  strain  invariants 
computed  by  explicitly  modeling  the  fiber  and  matrix  materials. 
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Failure  criteria  at  mid-ply 


Normalized  distance  to  free  edge 

Figure  1:  Plot  of  failure  criteria  versus  distance  to  free  edge  at  mid-ply  of 

both  the  0  and  90  degree  plies 

The  periodic  homogenization/de-homogenization  methods  are  upper  bounds  on 
the  actual  strain  in  the  explicitly  modeled  fiber/matrix  composite,  illustrating  that  the 
multi-scale  methods  give  accurate  results  when  the  microstructure  is  uniformly  periodic. 

As  one  approaches  the  0°/90°  interface,  these  methods  based  strictly  on  a 
uniformly  periodic  microstructure  are  no  longer  able  to  capture  the  trends  of  the 
explicitly  modeled  solution.  Figure  2  and  Figure  3  show  the  strain  invariant  plots  0.5 
fiber  diameters  above  the  ply  interface.  Note  in  Figure  3  how  the  second  strain  invariant 
in  the  full  numerical  simulation  exceeds  the  bound  from  the  multi-scale  bound  prediction 
assuming  uniformly  periodic  microstructure.  It  is  not  surprising  that  problems  would 
arise  in  this  region  since  the  microstructure  shifts  abruptly  from  a  0°  fiber  geometry  to  a 
90°  fiber  geometry. 

To  fix  ideas  for  the  asymptotic  theory  near  a  ply  interface,  we  illustrate  the 
approach  for  a  pre-stressed  laminate  made  up  of  several  plies  with  each  ply  containing  a 
different  periodic  array  of  fibers.  We  focus  on  a  domain  of  interest  "S”  inside  the 
composite.  The  characteristic  length  of  the  fiber  size  and  spacing  relative  to  the  length 
scale  of  the  laminate  is  given  by  s.  A  domain  of  interest  containing  the  interface  between 

two  plies  is  shown  in  Figure  4.  The  elastic  displacement  and  strain  inside  the  composite 

8  8 

laminate  are  written  as  u  and  e  respectively.  Next  we  replace  the  local  elastic  properties 
inside  each  ply  with  the  effective  elastic  tensor  associated  with  each  ply.  We  then 

H 

consider  the  elastic  strain  field  e  (x)  associated  with  the  laminate  composed  of  plies 
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modeled  by  their  effective  elastic  stiffness  properties.  The  new  high  fidelity  asymptotic 
representation  for  the  actual  strain  inside  the  domain  of  interest  S  is  given  by 

e*  (x)  =  Pe  (x)eH  (x)  +  e{rj e  (x))+)(s) ,  (1) 

8  8 

where  the  tensor  valued  function  P  and  vector  field  r\  are  determined  by  local  boundary 

8 

value  problems  formulated  over  the  domain  of  interest  S.  The  tensor  P  is  computed  in 

e,ij 

terms  of  the  local  elastic  fluctuation  w  that  is  periodic  on  S  and  solves 

div(C£  ( x)(eeJj  (x)  +  eu)  =  div{CE  (x)e,J ).  (2) 

E 

Here,  the  effective  elasticity  tensor  for  each  ply  is  given  by  C  (x),  the  local  elastic 

£  ij  £,ij 

property  for  the  fiber  and  matrix  phase  is  given  by  C  (x),  e  =8  5  and  e  is  the  local 

£,ij 

strain  inside  S  is  associated  with  w  and 

Pkf(x)  =  eskf(x)  +  el  (3) 

8 

The  field  r\  is  periodic  on  S  and  captures  the  local  strain  fluctuation  due  to  pre¬ 
stress.  The  local  stress  free  strain  takes  different  values  inside  fiber  and  matrix  phases  and 

8  8 

is  denoted  by  e  .  The  function  r\  is  the  solution  of 

div(Ce  (x)(e(r]e)  +  es  (x))  =  div(HE(x )) ,  (4) 

E  £ 

where  the  effective  thermal  expansion  coefficient  for  each  ply  is  given  by  H  (x)  and  e(r| ) 

8 

is  the  strain  associated  with  r|  .  Here  it  is  noted  that  all  differential  equations  are 
interpreted  in  the  weak  sense. 

Figure  5  illustrates  the  results  of  the  new  asymptotics  applied  to  the  previous 
0/90/0  free  edge  problem.  The  smooth,  black  lines  indicate  the  full  numerical  simulation 
explicitly  modeling  the  fiber  an  matrix  phases,  the  dashed,  blue  lines  portray  the  bound 
from  the  uniformly  periodic  microstructure  asymptotics,  and  the  dotted  red  line  displays 
the  bound  using  the  new  interface  asymptotics.  It  is  clear  that  the  previous  bound  is 
exceeded  by  the  full  simulation,  while  the  new  interface  bound  models  the  behavior  quite 
well.  This  work  was  submitted  for  publication  to  the  SIAM  Journal  of  Multiscale 
Modeling  &  Simulation  in  October,  2006. 
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Figure  2:  Plot  of  dilatational  strain  invariant  0.5  fiber  diameters  above  ply  interface 


Number  of  periods  to  free  edge 

Figure  3:  Plot  of  deviatoric  strain  invariant  0.5  fiber  diameters  above  ply  interface 
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Figure  4:  Domain  of  interest  containing  the  interface  between  two  plies  with  fiber 

orientation  at  0  and  90  degrees. 
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Number  of  periods  to  free  edge 


Figure  5:  Deviatoric  strain  invariant  at  ply  interface 
Task  2:  Comparison  of  failure  criteria  for  the  fiber/matrix  level 

This  task  was  not  originally  planned  to  be  a  part  of  the  work.  After  the  interface 
issue  was  addressed,  the  focus  was  to  optimize  the  fiber  orientation  angles  for  a  scarf 
repair  patch.  Before  we  proceeded  with  the  optimization  routine,  however,  we  had  to 
pick  a  failure  criteria  to  use  on  the  micro  level  (level  of  the  fiber  and  matrix).  We  used 
two  micro  level  failure  criteria  in  this  study.  The  first  was  stress  invariants  in  which  the 
dilatational  and  deviatoric  critical  values  were  deemed  independent  of  each  other.  In 
other  words,  no  interaction  was  assumed  between  the  two  failure  modes.  Second,  we 
used  the  criteria  for  metal  yielding  by  Doyoyo  and  Wierzbicki  [3].  We  note  that  the 
Raghava  criterion  ([4])  is  exactly  equal  to  the  Doyoyo-Wierzbicki  criterion  with  the 
shape  parameter  r\  set  equal  to  one. 

Failure  envelopes  were  computed  for  two  different  material  systems  under 
different  loading  conditions.  The  first  comparison  between  the  failure  criteria  was  an  in- 
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plane  shear-transverse  tension  test  on  E-Glass/LY556  composite.  For  calibration  of  the 
stress  invariant  criteria  (S-I),  we  used  90°  transverse  tension  to  back  out  the  critical 
dilatational  invariant  value.  The  critical  deviatoric  invariant  was  determined  by  pure 
shear  on  a  90°  specimen.  For  the  Doyoyo-Wierzbicki  criterion  (D-W),  we  calibrated  by 
pure  tension  and  pure  compression  on  the  90°  specimen,  then  adjusted  the  r|  parameter  to 
match  the  shear  point.  Figure  6  displays  the  experimental  and  predicted  failure  envelopes 
of  the  E-glass/FY556  composite.  Experimental  data  was  taken  from  [4,5].  Predictions 
were  done  with  and  without  prestress.  Because  of  the  low  stress-free  temperature  of  the 
composite,  little  difference  is  seen  in  the  analytical  predictions  with  and  without 
prestress.  It  is  clear,  however,  that  the  independence  of  the  failure  in  the  stress  invariant 
criteria  leads  to  a  poor  prediction  in  the  shear-tension  regime.  The  DW  criterion  appears 
to  capture  the  trend  much  more  accurately. 


Figure  6:  Tension/shear  failure  envelope  of  E-glass/LY556  composite 

The  second  failure  predictions  were  made  for  AS4/3502  angle-ply  composites 
loaded  under  compression.  Experimental  data  was  obtained  from  [6].  The  critical 
dilatational  invariant  value  for  the  stress  invariant  criteria  was  obtained  from  a  transverse 
tension  test.  The  deviatoric  critical  value  was  taken  from  the  ±45  angle  ply  under 
compression.  Calibration  for  the  Doyoyo-Wierzbicki  criterion  was  accomplished  using 
the  90  degree  specimen  under  tension  and  compression.  The  r\  parameter  was  then 
adjusted  to  match  the  ±45  angle  ply  experimental  data.  Once  again,  predictions  were 
made  with  and  without  prestress.  The  AS4/3502  composite  has  a  larger  stress-free 
temperature  than  the  E-glass/LY556  composite,  thus  differences  between  the  predictions 
are  not  negligible  in  this  case. 

Figure  7  displays  the  experimental  results  and  analytical  predictions  of  the  angle 
ply  data.  It  is  clear  that  the  D-W  criterion  performs  better  than  the  stress  invariant 
criteria.  It  should  be  noted,  however,  that  neither  criteria  accurately  capture  the  trend 
once  the  ply  angles  are  smaller  than  30  degrees.  It  is  believed  that  inside  30  degrees, 
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other  failure  modes  (such  as  fiber  buckling)  are  affecting  the  laminate  failure.  The  D-W 
criterion  with  prestress  shows  a  marked  improvement  over  the  D-W  prediction  without 
prestress  on  the  ±60  degree  specimen. 


Angle  Ply  (+/-  0) 


Figure  7:  Failure  of  AS4/3502  angle  ply  laminates 

Conclusion 

We  note  that,  at  present,  we  have  not  optimized  the  ply  angles  for  the  scarf  joint. 
The  micromechanical  failure  criteria  study  points  us  toward  interactive  criteria  for 
predicting  failure  in  the  matrix  phase  of  the  composite.  Although  the  Doyoyo-Wierzbicki 
criterion  appears  promising  as  a  micro-level  failure  criterion,  more  validation  is  needed  to 
test  the  situational  robustness  of  the  criterion. 

We  conclude  by  pointing  out  that  the  convergence  implied  by  (1)  holds  for  points 
inside  the  prescribed  domain  of  interest  S.  The  new  asymptotic  theory  shows  that  the 
approximation  improves  when  the  length  scale  of  the  heterogeneity  is  small  relative  to 
the  domain  of  interest.  Heuristically  this  implies  that  the  approximation  improves  if  the 
size  of  the  domain  of  interest  is  ' '  large  enough’  ’  with  respect  to  the  heterogeneity.  One 
can  think  of  the  domain  of  interest  as  being  a  particular  choice  of  Representative  Volume 
Element  (RVE).  The  question  of  selecting  the  proper  RVE  size  for  deterministic  and 
random  media  is  an  active  area  of  investigation  [7,8]  and  has  direct  impact  on  the  choice 
of  local  enhancement  of  the  FEM  seen  in  multi-scale  numerical  methods  [9-17].  The 
framework  given  by  the  asymptotic  expansion  (1)  provides  a  new  mathematical  context 
for  the  investigation  of  the  effect  of  the  location  and  size  of  the  RVE  on  the  fidelity  of  the 
approximation  and  choice  of  FEM  enrichment  for  multi-scale  numerical  methods. 
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